Stretched Exponential Relaxation on the Hypercube and the Glass Transition 
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We study random walks on the dilute hypercube using an exact enumeration Master equation 
technique, which is much more efficient than Monte Carlo methods for this problem. For each 
dilution p the form of the relaxation of the memory function q(t) can be accurately parametrized 
by a stretched exponential q(t) = exp( — (t/r)' 3 ) over several orders of magnitude in q(t). As the 
critical dilution for percolation p c is approached, the time constant r(p) tends to diverge and the 
stretching exponent /3(p) drops towards 1/3. As the same pattern of relaxation is observed in wide 
class of glass formers, the fractal like morphology of the giant cluster in the dilute hypercube is a 
good representation of the coarse grained phase space in these systems. For these glass formers the 
glass transition can be pictured as a percolation transition in phase space. 
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Complex systems generally show strongly non- 
exponential dynamics. In 1854 R.Kohlrausch used a 
phenomenological expression q(t) = Cexp(— (t/r)") to 
parametrize polarization decay data in Leiden jars |Q. 
Rediscovered more than a century later, again in the 
context of dielectric relaxation 0, this "stretched ex- 
ponential" or KWW expression has become ubiquitous 
in phenomenological analyses of relaxation data, experi- 
mental or numerical || . Analytical arguments have been 
given as to why certain hierarchical or trap model sys- 
tems show KWW relaxation [Bj^j, but there has always 
been a school of thought which considers that in the con- 
text of real glasses this expression is nothing more than 
a convenient fitting function of no fundamental signifi- 
cance. 

An alternative mechanism for KWW relaxation can be 
provided by a complex closed space approach. For ran- 
dom walks on a high dimensional critical percolation clus- 



ter in Euclidean space |]6[, it is known that (r 
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It was conjectured |7[] that random walks on a dilute 
hypercube in high dimension would necessarily lead to 
stretched exponential relaxation with a limiting value of 
/3 = 1/3 at percolation. Here we use an exact enumer- 
ation Master equation method which provides results of 
high precision for this problem. We find that the KWW 
functional form accurately fits data extending over a very 
wide range of q(t) values (from about 0.5 to 10 -5 ), and 
that j3 tends to 1/3 at percolation, as predicted. These 
results demonstrate that in the context of glassy systems 
the stretched exponential is not just an empirical fit- 
ting function, but that it has a precise mathematical and 
physical significance. A predictive generic picture for the 
glass transition can be obtained by identifying the glass 
transition with a percolation transition in phase space. 

Imagine a hypercube in high dimension N with a frac- 
tion p of its sites occupied at random. Clusters are de- 
fined as sets of occupied sites having one or more occu- 
pied sites as neighbors. It has been proved mathemati- 
cally [^) that there is a critical "percolation" concentra- 
tion p c given by 
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Where a = 1/(N — 1). For p > p c there exists a giant 
spanning cluster while for p < p c there exist only small 
clusters with less than TV elements. 

Now consider the relaxation due to random walks on 
the giant cluster of sites, a strictly mathematical problem 
which apparently has not been solved analytically. For 
a given realisation of the partially occupied hypercube 
with p > p c we can define a random walk among sites on 
the giant cluster. The walker starts at any such site i . 
A site j near neighbor to i is drawn at random. If j is 
on the giant cluster (and so "allowed" ) the walker moves 
to j. Otherwise the draw is repeated until an allowed site 



is found. Each draw, successful or not, is considered one 
time step. The procedure is iterated. 

We identify the distance Ha- between sites i and k on 
the hypercube with the Hamming distance, which is just 
the minimal number of moves needed to go from i to k on 
the full hypercube. The value of the normalized memory 
function q n (t) after time t for a given walk n starting 
from i Q and arriving at k n after time t can be defined by 
(2H ikri(t) - N)/(2N). The definition is identical to that 
of the autocorrelation function relaxation for the N Ising 
spins. The value averaged over many walks will go to 
zero at long t. 

Relaxation in the dilute hypercube has already been 
studied numerically by Monte Carlo techniques [p|,^o|. 
In the brute force Monte Carlo approach taking a mean 
over independent walks, the statistical noise becomes im- 
portant at long t, limiting accuracy flO]] . The exact enu- 
meration considers a Master equation to study the time 
evolution of the entire probability distribution for the 
walker after t steps, p(t), which we will call the state 
vector. Each vertex of the hypercube is associated to an 
integer < i < 2 N — 1. At t = the walker is localized 
on a single summit i a on the hypercube; the probability 
distribution then diffuses over the system at each time 
step following the equation: 



Pi(t)=pi(t-1) 



(2) 



- l)W(j - t) - Pj (t - l)W(i -> j) 



where W(i — > j) represents the transition probability 
that is given by: 

,. ( if i and j are allowed first neighbours 

^ ~~ * ' 1 otherwise 



Equation g can be rephrased as: 

pit) = Fp(t - 1) 



(3) 



(4) 



where F is the linear evolution operator. Our numerical 
algorithm catalogues all sites on the giant cluster for each 
particular realization of the hypercube, and then equa- 
tion (|^) is iterated for one particular starting point. Close 
to p c where time scales are long and there are fewer sites, 
it is more efficient to diagonalize the evolution operator 
F. 

By explicitly solving the master equation we obtain ex- 
act results (to within numerical rounding errors) for each 
combination of one realization of the hypercube occupa- 
tion, and one given starting point on the giant cluster. 
There is no statistical "noise" , and by averaging over a 
moderate number of independent samples a mean q(t) 
curve can be obtained, lying very close to the infinite 
ensemble average even to long times. 
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In practice calculations were done on dimension TV = 
16 for values of p from 0.5 to 0.073 (which is close to p c ). 
At least 100 samples were used at each p. N must be 
large; the present value was limited by computer memory 
considerations. 

We expect three relaxation regimes a priori. At short 
times q(t) will behave as 1 — at where a is the probability 
that a step will be made at a given attempt. Short time 
decay will thus be exp(— at). For very long t finite size 
effects will set in (the number of sites is finite for finite 
N) and a second crossover back to exponential relaxation 
will occur. Between these two limiting regimes, as the 
system explores the labyrinthine geometry of the giant 
cluster there will be the slow relaxation regime which 
interests us. 

The numerical data obtained together with the 
stretched exponential fits are shown as log q(t) against 
log(i) in Figure (|l|). The normalization parameter C is 
always close to 1. It can be seen immediately that the 
fits are of excellent quality and that they extend down to 
q(t) values as low as 10~ 5 . 

There are various methods of exhibiting this sort of 
data in order to make stringent tests of the functional 
form of q(t). For example, Figure @ shows log[— logg(i)] 
against log(t). In this plot perfect pure or stretched ex- 
ponentials should be straight lines at long t, with the 
pure exponential having slope one. Over a very wide in- 
termediate time regime, at each p the functional form 
of q(t) is indistinguishable from a stretched exponential. 
As p approaches p c the time scale r tends to diverge and 
the stretching exponent tends to near 1/3. The results 
strongly confirm the earlier conjectures and numerical 
work §g|g. 

In all relaxation models including the present one the 
shape of the decay can be related formally to a particular 
distribution of relaxation times of independent modes of 
the system. In trap models which have been studied in- 
tensively |p|,^|Jl2| individual non-interacting walkers fall 
into random traps, giving stretched exponential decay of 
the number of surviving walkers in the appropriate lim- 
its. The present model is not in the class of models with 
spins relaxing independently at different rates, and the 
initial site is not a privileged configuration. All the relax- 
ations of the single spins are coupled together implicitly 
through the complex morphology of the giant cluster. As 
the spins are strongly interacting, it is essential not to 
confuse the modes with the individual elements (spins) 
which are relaxing. Thus for the present calculations the 
effective number of spins is small (N = 16) but the num- 
ber of independent modes is much bigger: it is equal to 
the number of eigenstates of F, i.e. to p2 N , typically 
of the order of 10 4 modes for N — 16. The mode spec- 
trum is discrete and by definition has upper and lower 
limits. For times less than the minimum characteristic 
time, relaxation will be exponential, and for times much 
longer than the (N and p dependent) maximum charac- 



teristic time there must again be a second exponential 
regime, the finite size limit discussed above. The short 
time regime can be seen on all the curves, for q(t) values 
above about 0.5; the begining of the long time tendency 
to exponential decay is only visible for the lowest values 
of p where the number of modes is smaller and where the 
calculations have been taken to very long t. If calcula- 
tions could be carried out for much larger N the ultimate 
exponential regime would only appear at extremely long 
times and small q(t). 

In a sense the present model expresses concretely the 
physical picture proposed by Palmer et al [Q where the 
relaxation of each element depends on its instantaneous 
environment, but in contrast to the mode relaxation 
time distribution is not injected "by hand" but emerges 
naturally as a consequence of the fractal-like closed space 
morphology of the giant cluster, with no adjustable pa- 
rameters of any kind. It is important that the present 
approach not only leads necessarily to the stretched ex- 
ponential functional form, but provides an explicit quan- 
titative relation between the time scale and the stretch- 
ing. As p drops towards p c the time scale r(p) gets pro- 
gressively longer (a divergence at p c in the very large N 
limit). Concomitantly 0{jp) decreases from 1 at large p 
towards a limit of 1/3 at p c , (figure (||)). Both effects 
reflect the increasing sparseness and complexity of the 
giant cluster with decreasing p. The limiting value of 
1/3 for when r diverges is a consequence of the "per- 
colation fractal" morphology of the sparse giant cluster 

We can now compare heuristically the model data with 
examples of numerical and experimental relaxation re- 
sults. The autocorrelation function relaxation in the 3d 
bimodal Ising spin glasses has been studied numerically 
to high precision [g_3[. The long time relaxation func- 
tion above the ordering temperature is of stretched ex- 
ponential form with an exponent which tends to 1/3 
within numerical accuracy at the temperature at which 
the times scale r diverges. Relaxation in other spin 
glasses are of the same limiting form, independently of 
space dimension or of the type of spin-spin interaction 
|]l4| . This pattern of behaviour is not restricted to spin 
glasses. For instance the relaxation of a colloid glass 
former Ea] , a system having an entirely different micro- 
scopic mechanism for glassyness, again shows stretched 
exponential decay with tending to 1/3 as r diverges 
with concentration because of steric hindrance. A large 
number of polymer glass formers also show a similar char- 
acteristic relaxation pattern [jfq]. 

What is the logical connection between the dilute hy- 
percube random walk and relaxation in glasses ? The 
ensemble of all possible configurations of a thermody- 
namic system form a high dimensional closed space. The 
available phase space at temperature T can be consid- 
ered at the microcanonical level as the set of configura- 
tions having the appropriate energy for that temperature, 
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E(T). This available phase space becomes sparser as T 
decreases. Phase transitions correspond to discrete qual- 
itative changes in the morphology of the available phase 
space with temperature; thus at a standard ferromag- 
netic second order transition T c the phase space splits 
into two. Also quite generally, relaxation is just the con- 
sequence of the random walk of the configuration point 
of the whole system in the available phase space, and its 
form is necessarily a reflection of the morphology of this 
phase space. Explicitly the N dimension hypercube is ex- 
actly the total phase space of an N spin Ising model; the 
spin by spin relaxation of a coupled N spin Ising system 
can be mapped directly onto a random walk of the config- 
uration point on the thermodynamically available sites of 
the N dimension hypercube |]l3| . For systems with more 
complicated total phase spaces than the hypercube, the 
same argument applies mutatis mutandi. 

A non-exponential relaxation means a complex phase 
space; the striking ressemblence between the dilute hy- 
percube relaxation pattern and the relaxation observed 
in the numerical spin glasses or experimental glasses we 
have cited necessarily implies similarities at the level of 
phase space geometry. Thus, at short times (fine grained 
phase space) the relaxation will depend on the details 
of the physics of each system, but at moderate and 
long times (coarse grained phase space) these systems 
all appear to have the same specific percolation-fractal- 
like phase space morphology with its characteristic relax- 
ation signature, the precursor of a phase space percola- 
tion breakdown at the glass transition. By continuity, be- 
low the percolation transition the phase space would split 
into many inequivalent clusters, meaning for the thermo- 
dynamic systems a transition to a frozen glassy state. 
The glass transition corresponding to this description is 
a continuous transition, but in an entirely different class 
from standard second order transitions. 
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FIGURE CAPTIONS 
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FIG. 1. Decay of the autocorrelation function q(t) on a 
log-log plot for different values of p, as listed in the inset. The 
solid lines correspond to stretched exponential fits, with f5(p) 
and r(p) as indicated in Figure 3. The error bars correspond 
to an estimate of the uncertainty of the points due to limited 
sampling. 
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FIG. 2. A more stringent test of the stretched exponential 
behavior of q(t) is the log [-\og(q(t)) ] against log t plot. The 
different values of p are listed in the inset. The solid lines 
correspond to a stretched exponential fit, with (5 and r as 
indicated in Figure 3 
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FIG. 4. Experimental data from light scattering mea- 
surements on a polystyrene colloid, taken from Bartsch et 
al ref [15]. In the publication the relaxation curves were 
parametrized using the KWW form. The relaxation was mea- 
sured as a function of the colloid volume fraction <f>. The 
critical value <f> c where a gel forms is about 0.69. 




FIG. 3. Relaxation time r(p) (circles) and stretched ex- 
ponential exponent f3(p) (squares) against p. As p — * p c , r 
diverges while (3(p) approaches 1/3. 
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